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The scattering of small bodies in planetary systems: 
constraints on the possible orbits of cometary material. 
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ABSTRACT 

The scattering of small bodies by planets is an important dynamical process in 
planetary systems. In this paper we present an analytical model to describe this process 
using the simplifying assumption that each particle's dynamics is dominated by a 
single planet at a time. As such the scattering process can be considered as a series 
of three body problems during each of which the Tisserand parameter with respect 
to the relevant planet is conserved. This constrains the orbital parameter space into 
which a particle can be scattered. Such arguments have previously been applied to the 
process by which comets are scattered to the inner Solar System from the Kuiper belt. 
Our analysis generalises this for an arbitrary planetary system. For particles scattered 
from an outer belt directly along a chain of planets, based on the initial value of the 
Tisserand parameter, we find that it is possible to (i) determine which planets can 
eject the particles from the system, (ii) define a minimum stellar distance to which 
particles can be scattered, and (iii) constrain range of particle inclinations (and hence 
the disc height) at different distances. Applying this to the Solar System, we determine 
that the planets are close to optimally separated for scattering particles between them. 
Concerning warm dust found around stars that also have Kuiper belt analogues, we 
show that, if there is to be a dynamical link between the outer and inner regions, then 
certain architectures for the intervening planetary system are incapable of producing 
the observations. We speculate that the diversity in observed levels of warm dust may 
reflect the diversity of planetary system architectures. Furthermore we show that for 
certain planetary systems, comets can be scattered from an outer belt, or with fewer 
constraints, from an Oort cloud analogue, onto star-grazing orbits, in support of a 
planetary origin to the metal pollution and dustiness of some nearby white dwarfs. 
In order to make more concrete conclusions regarding scattering processes in such 
systems, it is necessary to consider not only the orbits available to scattered particles, 
but the probability that such particles are scattered onto the different possible orbits. 



1 INTRODUCTION 



The scattering of small bodies is an important dynamical 
process in many planetary systems. One classic example is 
the population of small bodies close to the Sun, many of 
which originate further out in the Solar System, from where 
they were scattered inwards. Near-Earth asteroids (NEAs) 
originate in the asteroid belt. Many left the belt after be- 
ing destabilised by resonances with Jupiter and then scat- 
tered by the terrestial planets (|Morbidelli et al]|2002l '). Vis- 
ible comets are either objects s cattered inwards from th e 
Kuiper belt or the Oort cloud (|Levison fc DuncanI 119971 ). 
The scattering of small bodies has not been considered in 
detail for extra-solar planetary systems, mainly due to the 
lack of constraints on the structure of the planetary system. 
There is, however, evidence for small bodies in many extra- 
solar planetary systems. Dust belts, known as debr is discs 
are s een around hundreds of main sequence stars (jWyattl 
I2OO8I ). Observations, particularly resolved images, suggest 



that debris discs interact with planets jKalas et al.l l2005l : 
iMoerchen et al.ll201ll :[G reaves et al.ll2005r ). etc. Assuming a 
similar nature to our Solar System, it is reasonable to ex- 
pect that scattering in these systems can also result in a 
comet-like population. The expected level and distribution 
of this comet population may differ substantially from the 
Solar System, depending on the individual planetary system 
architecture. 

Evidence of such a comet-like population may exist 
from observations of wa rm dust discs aro u nd a handfu l 
of main sequence stars (IWvatt et al.l I2OO5I : iGaidod [1999I : 
iBeichman et al.ll2005l : [Song et al.ll2005l ). Comets or asteroids 
in the position of the observed dust belts have a short life- 
time against collisions and drag forces. They cannot have 
existed for the entire main sequen ce lifetime in their ob- 
served position (|Wvatt et al.|[2007l ). One possible explana- 
tion is that the material originated in a cold, outer belt. 
It could be that we are observing a comet-like population. 
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that is continuously replenished from scatter ing of mate- 
rial fr om the outer belt by intervening planets (|Wvatt et alj 
[20o3)- Alternatively, it could be a transient ever it, resulting 
from the stochastic collision of two larger bodies (|Song et alJ 
[iooi), maybe in a similar manner to the impact that formed 
the Earth-Moon system. Or, material could be transported 
inwards from the outer belt during a LHB type event 
llBooth et al.ll2009l ) or by drag forces (|Reidemeister et alJ 
120111). 

Another piece of evidence for the scattering of mate- 
rial in exo-planetary systems comes from observations of 
evolved stars. 2 5% of DA white dwar fs show unexpected 
metal pollution jZuckerman et al]|2003h . whilst 1-3% of DA 
white dwarfs with cooling ages less than O.SGyr have excess 
emission in the inf ra-red consistent with a close-in dust disc 
l|Farihi et al.ll200^ ). The comp osition of the pollu ting mate- 
rial closely resembles planets (|Klein et al.ll201oh and there 
is good evidence tha t it is not materia l accreted from the 
inter-stellar medium (|Farihi et al ] |2010al ). The best theories, 
therefor e, suggest that it originates in an outer planetary 
vstem (I Jurall 2003': 'Gansickc ct al.' '200^; 'Kili c et all l2006l: 



von Hippel et al. 2007 : Farihi ct al.,,2009, . 2010b|; iMelis et all 
As the star loses mass on the giant branch, dynami- 



cal instabilities can be indu ced in the outer planetary system 
l|Debes fc Sigurdsson|[20o3 ). These can lead to comets or as- 
teroids being scattered by interior planets onto star-grazing 
orbits, where they are tidally disrupted. Material from the 
tidally disrupted asteroids or comets forms the observed 
discs and accretes onto the star. The ability of evolved plane- 
tary systems to scatter comets or asteroids onto star-grazing 
orbits requires further detailed investigation, although previ- 
ous work ha s considered an Port cloud origin of the scatter ed 
bodies e.g. (|Alcock et al.ll 19861 : lOebes fc Sigurdssonll2002l ). 

In this work the scattering of small bodies in 
an arbitrary planetary system is investigated. N-body 
simulati ons are typically used to model such scat- 
tering jLevison fc Duncanl Il997l : iHorner fc Joned l2009l : 
iHolman fc WisdomI 11993^ A deeper understanding of the 
general properties of such scattering can, however, be 
achieved using analytical arguments. Simulations of scat- 
tered Kuiper belt objects have found that the scattering pro- 
cess can be approximated as a series of three-body problems, 
as th e scattered bodies are p assed from one planet to the 
next iLevison fc Duncanlll997l ') . While such particles are un- 
der the influence of one of the planets, their dynamical evo- 
lution can be approximated by the circular restricted three- 
body problem in which the orbits of the p articles must be 
such that their Tisserand parameters, Tp, (|Tisserandl 1 18961 : 
iMurrav fc Dermott|[l999l ') are conserved, where 



T, = ^+2,/(l^£>cos(/), 



(1) 



where a, e, I are the comet's semi-major axis, eccentricity 
and inclination and is the planet's semi-major axis. This 
conservation is so fundamental to come tary dynamics that 
it is used to classify cometary orbits (| Horner et al.l l2003l : 
iGladman et~al]|2008l ). 

In this work we use the conservation of the Tisserand 
parameter to constrain the orbits of scattered particles in 
a planetary system with an arbitrary configuration. In sec- 
tion [2] we discuss how planetesimals are scattered from an 
outer belt, in an otherwise stable planetary system. We then 



outline our constraints on the orbits of particles scattered by 
a single planet in section |3l which we extend to two planets 
in section U and arbitrarily many planets in section [S] In 
section [6] we consider the application of this analysis to our 
Solar System, systems with warm dust discs and polluted 
white dwarfs. 



2 SCATTERING OF PLANETESIMALS 

During the planet formation process, a planet that forms in 
a disc of planetesimals, will swiftly clear a zone around it, 
both by scattering processes and resonant interactions with 
the planet. Analytically the size of the planet's cleared zone 
can be approximated. Criterion for the overlap of mean mo- 
tion resonances determine a re gion around the planet within 
which orbit's are chaotic (|Wisdom,1980l ). whilst the Jacobi 
constant can be used to det ermine the zone within wh ich or- 
bits can be planet crossing jCladman fc Duncanlll990l ). Sim- 
ulations have shown that Neptune clea r s such a zone in less 
than 10^ yr (|Levison fc Duncanl Il993l : iHolman fc WisdomI 
Il993h . but more generally one might expect 1,000 conjunc - 
tions for this clearing to take effect (jDuncan et al.lll989l ). 
Material removed from this region may be ejected, whilst 
some fraction remains on bound, eccentric orbits, with peri- 
centres close to the planet's orbit, forming an analogue to 
Neptune's scattered disc. After many scatterings some of this 
material may reac h far enough from the star to interact with 
the Galactic tide (|Tremain 3 119931) and eventually populate 
an analogue of the Oort cloud. 

Planetesimals outside of this zone could in principle 
be long term stable. However, N-body simulations of Nep- 
tune and the Kuiper belt find that Kuipe r belt objects 
are s ti ll scattered by Neptune a t late times (iDuncan et ahl 
19951: IHolman fc WisdonJ Il993l: iLevison fc Duncanl 119971 

m 



Emel'vanenko et al.ll2004l : fMorbideUilll997f ). The Kuiper beh 
has a complicated structure of stable and unstable regions. 
The gravitational effects of Neptune and the inner planets 
result in the overlap of secular or mean-motion resonances 
produci ng thin chaotic regions, within the otherwise s table 
region (|Kuchner et al.l |2002| : iLvkawka fc Mukail [iooi) and 
small unsta ble regions within otherwise stable mean motion 
resonances (iMoons fc MorbideUi|[l995l : iMorbideUi fc Moond 
1 19951 : iMorbidellil Il997l ). Objects may diffuse cha, otically 
from stable to unstable regions (jMorbidellil |2005| ). This 
has been shown to occur for Neptu ne's 3:2 



process 



and 2:1 r esonances, amon gst others (jM orbidellil 



Tiscareno fc Malhotra 2009; iNesvornv fc Ro ig 200Q,, 



1997 



2001 



de Eha et al.l |2008| ). Objects leaving mean motion reso- 



nances in the Kuiper belt, in this way, may be the main 
source of Neptu ne encountering ob jects at the age of 
the solar system (|Duncan et al.l 19951 ). Many of these ob- 
jects are scattered into the inner planetary system, and 
could be the source o f Centaurs or Ju p iter Family comets 
ddi Sisto et all I2OIOI: IMorbidellil 1 19971 : iLevison fc Duncanl 
ll997l : lHolman fc Wisdomlll993l ). 

The dynamical processes occurring in the Kuiper belt 
may well be applicable to exoplanetary systems with a sim- 
ilar structure, i.e. an outer planetesimal belt and interior 
planets. The outer be lt could be truncated by resonance 
overlap (| WisdomI 1 1980l ) . Most particles would then inhabit 
a predominately stable region exterior to this, containing 
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small regions that are unstable due to the overlap of secu- 
lar or mean motion resonances of the inner planets. Objects 
could diffuse chaotically on long timescales from the stable 
to unstable regions and be scattered by the outer planet. 
Some of these scattered objects could enter the inner plan- 
etary system, whilst some could be ejected. 

In our consideration of the dynamics of material scat- 
tered from the outer belt by interior planets, we find that 
this dynamics is strongly dependent on the initial value 
of the Tisserand parameter, with respect to the outermost 
planet. Therefore, it is important to consider the value of this 
parameter. Objects in the outer belt tend to have T > 3, 
whilst in order that an object be scattered by the planet, 
the Tisserand parameter must be less than 3. Hence, for 
particles that are scattered at late times due to chaotic dif- 
fusion into an unstable region, at the time of first scatter- 
ing, the Tisserand parameter would be expected to be close 
to 3. Simulations of our solar system found this to be the 
case (ILevison fc Duncan|[l997l ). Here, we consider the initial 
value of the Tisserand parameter of such particles as an un- 
known, with the expectation that objects scattered in the 
way described will have initial Tisserand parameter values 
just below 3. 

As a final note, we point out that not all objects scat- 
tered by Neptune originate from the cold Kuiper belt. The 
two other main sources are Neptune's scattered disc and the 
Oort cloud. It is possible that similar classes of objects ex- 
ist in exo-planetary systems, however, there is at present no 
evidence for exo-Oort clouds or scattered discs. The distribu- 
tion of the Tisserand parameter for such objects would differ 
significantly from those that leave the cold Kuiper belt, in 
particular for Oort cloud objects, where it is unconstrained 
and T < 2 is possible. Therefore, for clarity and simplic- 
ity, in this work we focus on the objects that originate in 
an outer belt and that are first scattered by the outermost 
planet, at the age of the system. 



3 SCATTERING BY A SINGLE PLANET 

Firstly we consider a system similar to that described in the 
previous section, with a single planet, labeled by subscript 
i on a circular orbit at at and an exterior planetesimal belt. 
We consider planetesimals scattered from the outer belt by 
the planet. We make the simplifying assumption that plan- 
etesimals only interact with the planet if their orbits directly 
cross the planet's orbit. This simplifies the following analysis 
and enables analytical limits to be easily derived. However, 
since in reality interactions will occur in a zone around the 
planet, care should be taken in rigorously applying any of 
the derived limits, in particular for more massive planets. 
This will be discussed further in 33 

3.1 Orbital constraints 

For a planetesimal with a given value of the Tisserand pa- 
rameter with respect to this planet, Ti, the potential orbits 
onto which it can be scattered are limited, no matter how 
many times it interacts with the planet. The Tisserand pa- 
rameter gives us no information about the probability for 
any given interaction to scatter a planetesimal onto a given 
orbit, nor the timescales for interactions to occur. It does. 



however, limit the orbital parameters of the planetesimals 
after the interaction, in terms of its pericentre, q, eccentric- 
ity, e and the inclination, /, of its orbit with respect to the 
planet's. These constraints can be represented by a 3D vol- 
ume in {q, e, I) space. A planetesimal, given an initial value 
of Ti, may not be scattered onto an orbit with parameters 
outside of this volume, in this simple example. 

This parameter space can be fully mapped out analyt- 
ically by re- writing Eq. [T] as 

T. = Ml^+2/(i±^cos(/), (2) 

and noting that if the planetesimal is to remain on a bound 
orbit, < e < 1, — 1 < cos(J) < 1 and q > must apply. 
In order that the particle is scattered by the planet its orbit 
must cross the planet's and thus Q > ai and q < ai must 
apply. Applying these constraints to Eq. (2] places analyti- 
cal bounds that define this 3D volume of permitted orbits. 
Given the difficulties in presenting a 3D volume, we instead 
present the 2D projection of this 3D volume onto the q — e 
plane, I — q plane and e — I plane, shown in Fig. [1] The 
analytical bounds are presented in Table. [T] 

3.2 Minimum pericentre 

Further examination of the q-e plot in Fig.[T]makes clear that 
planetesimals cannot be scattered further towards the star 
than a limiting value, qmin , determined by T. This value can 
be calculated using constraints on the orbital parameters, 
Q = ai and cos(/) = 1 (equivalent to the lower bound in the 
q-e plane). For 2 < T, < 3: 

qrmn _ -Tf + 2Ti + A - ^sfZ^^ 
a, - Ti-8 ■ 

qmin as a function of Ti is shown in Fig. (2] The eccentricity 
at qmin will be given by: 

eum = T-3 + 2^3~^. (4) 

For Ti < 2, the lines Q — ai and cos(/) = 1 (positive 
root) no longer cross and the parameter space in the q-e 
plane is no longer bounded by Q = Oi, rather by cos(/) = 1 
(both positive and negative root). Therefore qmin — >■ 0. This 
can be shown to be true by considering the derivatives of 
the lines: 

dq, d,q. , , 

-r\coBl=l,g^O > -r\Q = ai,g^O- (O ) 

ae ae 

Importantly this implies that the constraints on the 
pericentre that apply to the orbits of objects with T > 2 
are not applicable to those with T < 2; such objects can be 
scattered onto orbits with any pericentre. 

3.3 Ejection 

A single planet can also eject planetesimals, given a suit- 
able value of the Tisserand parameter. Unbound orbits {i.e. 
those with e > 1) are not included in the plots in Fig. [T] 
It is, however, possible to determine from the top panel of 
Fig. [1] those values of the Tisserand parameter for which 
the particles are constrained to bound orbits with e < 1. 
The most eccentric orbits are those with pericentre at the 
planet's orbit (g = ai), therefore substituting into Eq. [2j 
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Figure 1. The possible orbital parameters of particles scattered by a single planet, with a given value of the Tisserand parameter with 
respect to that planet, T^. This forms a 3-D parameter space, that is shown here projected onto the eccentricity-pericentre (g — e) plane, 
the inclination-pericentre {I — q) plane and the eccentricity-inclination (e — /) plane. The limits of the parameter space are defined 
analytically in the Table[T] The units of the pericentre is the planet's semi-major axis and the dotted black line in the top row of plots 
shows the line where the particle's semi-major axis is equal to the planet's (a = = a^). 



(qr = ai,e = 1,1 = 0°), we find that there is a limit 
on the Tisserand parameter such that only objects with 
can be ejected. This has been pre v iously calcu- 
lated in, amongst others, iLevison fc DuncanI (|l997f ). usme; 
the formulation for the Tisserand parameter of a parabolic 
orbit. It should, however, be noted that this only applies 
strictly for low mass planets. As the planet mass is increased, 
as does the zone of influence of the planet. If we assume 
that particles can interact with a planet if they are within 
a distance A from the planet, then a particle on a hyper- 
bolic orbit with Tp = 3 and g = |ap can still be ejected if 
g < flp-l- A. A will be a function of the planet's Hill's radius , 
Rh , for example if A ~ 2V3Rh (jdadman fc Duncanlll990l l 
Jupiter can eject particles with Tjup = 3, whilst Neptune 
cannot. 



4 SCATTERING BY TWO PLANETS 

Now consider a planetary system with an outer belt and two 
interior planets, both on circular orbits. Particles from the 
outer belt are scattered by the outer planet, 1. The main 
possible fates of such particles are ejection, collision with a 
planet or the star, further scattering interactions with this 
planet, or scattering by the inner planet, 2. Many scattered 
particles are scattered multiple times by the outer planet. 
It dominates their dynamics for a certain period of time, 
during which the Tisserand parameter, with respect to this 
planet, Ti, is conserved. At some point, the particle may be 
scattered onto an orbit that overlaps with the inner planet 
and it may be scattered by that planet. In such an inter- 
action the Tisserand parameter with respect to the inner 
planet, T2 would be conserved, rather than Ti. Depending 
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Table 1. The analytical boundaries on the parameter space constraining the potential orbital parameters of a particle scattered by a 
planet, where the initial value of the Tisserand parameter is Ti. All units are in terms of the planet's semi-major axis; = 1. For the 
cases where more than one limit is stated, the upper of the two applies. 
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Figure 2. The minimum pericentre for a test particle scattered 
by a single planet, as a function of the Tisserand parameter value, 
from Eq. |3] For Ti < 2, q^in ^ 0. 

on the new orbit, it is then hkely that the particle is re- 
scattered by the inner planet and for a certain period its 
dynamics will be dominated by that planet. 

We start by considering this simple situation where the 
particle is passed from planet 1 to planet 2. This is used 
to describe constraints on the orbits of scattered particles. 
We then consider the possibility that particles are scattered 
backwards and forwards between the two planets in sec- 
tion [331 

4.1 Orbital constraints 

For a particle scattered by the outer planet the Tisserand 
parameter, Ti, is conserved. The value of T\ constrains the 
orbits, (5, e, 7), of scattered particles to those shown in Fig.[T] 
that satisfy Eq. (2] Although only sets of the orbital param- 
eters, q,e,I, that satisfy Eq.[2]are allowed, the full range of 
possible values is given by: 

q e [gmm(Ti), 1] 

I £[Q,ImaATl)l (6) 

where qmin is given in Eq. [3l 



T;=2.9 
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Figure 3. The orbital parameter space as determined by the 
Tisserand parameter, in the eccentricity-periccntre plane, for Ti = 
2.9 (equivalent to the fourth plot on the top row of Fig. [TJ. The 
bounds on this space are between the dashed (cos 7 = 1) and 
dot-dashed lines {q = ai) and shown in green. The subset of this 
orbital parameter space that can interact with an inner planet 
placed at ai„ = O.Sa; is shown by the vertically hashed region, 
whilst the subset that could interact with an outer planet placed 
at aout = 1.5ai is shown by the horizontally hashed region. The 
dotted line shows a = ai, the solid line q = ai„ and the triple 
dotted dashed line, Q = aout- 

W = cos"'(Vri-2), (7) 
and 

e™^ = 3 - Ti + 2V3~^. (8) 

As mentioned earlier, if Ti > 2\/2, then Cmax > 1 and 
some orbits are unbound. 

The particle may interact many times with the outer 
planet, moving between orbits in this parameter set, until at 
some point it encounters the next planet, 2. Only a subset of 
the orbits specified by Ti can interact with the next planet, 
2. These are shown by the green filled area in Fig.[3]and are 
those orbits that cross the planet's, with g < 02 and: 

q G [qm.in{Ti),a2] 

e G [e,„t( — ), eiim(ri)] 
ai 
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Figure 4. The minimum periccntrc for a test particle scattered 
by two planets, as a function of the ratio of the inner planet's 
semi-major axis to the outer planet's semi-major axis. 



where, 



ai + 02 
ai — 02 



(9) 



(10) 



If < iUri ' then the set of orbital parameters with 
I = Imax do not cross the inner planet's orbit. This occurs if 
the second planet is inside the maximum in / as a function 
of q that occurs &t q = |^-^ai (see Fig. [1]). In which case I 
is constrained to be less than Unt rather than Imax, where 



lint = cos 



ri(i + ft) 



(11) 



Once the particle is scattered by planet 2, Ti is no longer 
conserved, instead the value of T2 when the particle is first 
scattered by planet 2 is conserved. The range of possible 
T2 values is determined by the initial value of Ti and the 
planets' orbits, specified by the ratio of the planet's semi- 
major axes, 1^. 

The minimum possible value that T2 can have occurs 
for particles on orbits with minimum pericentre {q = qmin), 
the correspond eccentricity (e = eum) and in the orbital 
plane of the planets (7 = 0°). It is given by: 



T2,: 



02(1 



+ 2 



(1 + enm)qrr 
02 



(12) 



where eu-m (Eq- SI and qmin (Eq. [Sjl are functions of Ti . 

Since the Tisserand parameter {T2) is a monotonically 
increasing function of q, T2 will be maximum for the orbit 
with the largest value of the pericentre, q, that still crosses 
the planet's orbit, i.e. q — 02. For the range of T2 values 
for orbits with g = 02, the minimum is at cos / = ±1 and 
e = ex, where from the top line of Table [1] 



ex 



i+2(^)^_(^)ri+2(^)«/%/ 
01 01 01 y 



2-K^)3_(^)Ti.(13) 

Ol Ol 



Hence, the maximum of T2 is given by: 
T2,max = (1 - ex) + 2^/1 + ex- 



For the next time period the dynamics of the particle 
is controlled by the second planet. It may be scattered once 
or many times. Yet again, the particle's orbit is constrained 
to orbital parameters, (g, e, J), specified by the value of T2 
and Eq. [2] This time, however, we consider the situation 
where only Ti and the planet's orbits are specified initially 
such that it is only known that T2 lies between T2.min and 
T2,,nax- The full range for the orbital parameters {q,e,I) is 
therefore specified by: 



i(r2,min) < g < 1 

^ e <^ emax(T2,niin^ 
^ / *C /maLC (/2,min ) ; 



(15) 
(16) 
(17) 



where qmin is given by Eq. [31 emax by Eq. [5] and Imax by 
Eq. [T] but as a function of T2,min rather than Ti. 



4.2 Constraints on which particles interact with 
the inner planet 

For specific planetary orbits, specified by the ratio of the 
planets' semi-major axes, and strict constraints on the 
initial value of the Tisserand parameter in the outer belt 
{i.e. Ti close to 3), the orbits of scattered particles may 
be constrained such that they never interact with the inner 
planet. This occurs when the minimum pericentre to which 
particles may be scattered by the outer planet is further 
from the star than the inner planet's orbit; gmm(Ti) > 02 
(Eq. |3]) or : 



-Tf -I- 2ri + 4 



- 4V3 - Ti 02 
8 '^01 



(18) 



(14) 



4.3 Minimum pericentre 

In Sec. 13.21 Eq. [S] we determined the minimum pericentre to 
which a single planet may scatter a particle. A similar calcu- 
lation may be made for two planets, assuming that particles 
are only passed once along the chain of planets. The min- 
imum pericentre will depend on the Tisserand parameter 
with respect to the outer planet, Ti and the ratio of the 
planets' semi-major axes, 

For a particle that is scattered by the outer planet, 
with a value of the Tisserand parameter with respect to 
that planet of Ti, if it is then scattered by the inner planet, 
the particle could have a range of possible values of the Tis- 
serand parameter with respect to the inner planet, between 
T2 {Eq.\Ui and T2,max (Eq.[l4|. S 1I1C6 Qmin (Eq. |3} is a 
monotonically increasing function of the Tisserand param- 
eter, the minimum pericentre for scattering by both plan- 
ets will be given by qmin{T2,min), where T2,min is the min- 
imum value of the Tisserand parameter (Eq. I12|) . If a par- 
ticle is to eventually be scattered inwards as far as possi- 
ble by the outer and inner planet, it must be passed from 
the outer to the inner planet with an orbit of eccentricity 
e = eumiT2,min) (Eq. [i]) and inclination, / = 0°. 

The minimum pericentre for a two planet system is 
shown in Fig. [4] as a function of the ratio of the planets' 
semi-major axes, This is calculated from Eq. [3] such 
that q = qmin{T2.min), where T2,min — T{q — ap,e = 
eiim{Ti),I = 0°), using Eq. [2] This has a clear minimum, 
which occurs at: 
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0>2,min 



{l + eHra{Tl)y^''qrmn{Ti) 

(l-ei™(Ti))2/3 



(19) 



where e;i„i and Qmin are the minimum pericentre and hmit- 
ing eccentricity for scattering by the outer planet, given by 
Eq. Eland Eq.E 

This means that the optimum configuration of two plan- 
ets in terms of their abihty to scatter particles as close to the 
star as possible, involves planets positioned in semi-major 
axis at a2,min and ai. It is interesting to note that the op- 
timum position for the inner planet is not as close to the 
star as the outer planet could possibly scatter particles i.e. 
?min(7i), but closer to the outer planet. This is because 
there is a balance between moving the inner planet closer to 
the star, decreasing 02, such that qmin is decreased directly 
or moving the planet further from the star, increasing 02, 
but decreasing T2 and thus qmin. Of course this does not 
include any information about the probability that the par- 
ticle is ejected or collides with the planet rather than being 
ejected. 



4.4 Further scattering 

Scattering is not confined to the forward direction. Particles 
may originate in the outer belt, be scattered inwards by the 
outer planet, passed onto the inner planet, and then scat- 
tered back outwards again to the outer planet. Constraints 
on which particles might re-interact with the outer planet 
can be determined using a similar procedure to that dis- 
cussed in the previous section (Sec. 14.1] ) for particles passed 
from the outer planet to the inner planet. 

The possible values for the orbital parameters of par- 
ticles scattered by the inner planet are determined by the 
value of the Tisserand parameter, T2. A subset of these or- 
bits cross the outer planet's orbit, those with apocentres 
outside of its orbit (Q > ai). For the example of an outer 
planet at ai — 1.5a2 and with T2 — 2.9, this subset is shown 
by the hashed region in Fig. (3] Each set of orbital param- 
eters in this region {q, e, I) will specify a possible value for 
the Tisserand parameter with respect to the outer planet, 
Ti. The minimum possible new value of Ti occurs at the 
maximum pericentre (q — 02), the maximum eccentricity 
{emaxiT^) Eq. [8]) and cos / — 1, such that: 



the particle is then scattered back outwards, the minimum 
value of Ti is 2.93 (Eg. I20|) . If the particle is then scattered 
back in, again from Fig. 3] this gives a new minimum peri- 
centre for scattering by the two planets of qmin ~ 0.12ai . Af- 
ter a further scattering backwards and forwards, qmin — 0; 
all constraints on the eccentricity and pericentre of the orbit 
are removed. Given sufficient repetitions this occurs for all 
pairs of planetary orbits, where the constraints on the Tis- 
serand parameter allow particles to be passed between them. 
Thus, the orbital parameter space available to scattered par- 
ticles can be greatly increased by repeatedly scattering them 
backwards and forwards. 

So far we have merely outlined the orbital parameter 
space available to particles and not discussed the probabil- 
ity for scattering particles into this space. This is in general 
beyond the scope of this paper, however, these have impor- 
tant implications for the passing of particles backwards and 
forwards between the two planets. Firstly, it is clear that the 
timescales for particles to be repeatedly scattered backwards 
and forwards between two planets will be long and therefore 
at any given time the probability will be higher that par- 
ticles have merely been scattered by the outer planet, or 
passed from the outer to the inner planet once. Secondly, al- 
though repeated passing of particles between planets greatly 
increases the range of orbital parameters available to such 
scattered particles, this does not mean that it is most prob- 
able for such particles to be scattered onto more extreme 
(higher eccentricity or inclination) orbits. In fact, if we were 
to assume that a particle has an equal probability of being 
scattered onto any of the orbital parameters available to it, 
it is most likely that the particle is scattered onto an orbit 
that retains a value of the Tisserand parameter close to its 
original value. It is only the few particles that are scattered 
onto extreme orbits, i.e. with low pericentre or high eccen- 
tricity/inclination, that have significantly reduced values of 
the Tisserand parameter when they are scattered by the 
next planet. Therefore, although it is possible that particles 
may be scattered onto extreme orbits, with low values of 
the Tisserand parameter, by being repeatedly passed back- 
wards and forwards between the planets, we anticipate that 
the probability for this to occur is low and we are therefore 
justified in focusing on particles scattered directly along a 
planetary system for the rest of the paper. 



' l,new ,min — 



ai(r2 



2^/3^^) 



+2 



a2 



{4-T2 + 2y/3~n)a2 



ffli 



(20) 



If there are a range of values for T2, the smallest {e.g. T2,min 
for Eq. I12p will give the lowest value of Ti^naw.min. The 
maximum value of T\ such that particles can still interact 
with the outer planet is 3, as for any scattering event. 

If the particle is scattered backwards and forwards mul- 
tiple times this procedure may be repeated to determine 
the full range of Tisserand parameter values and potential 
orbits. T'i,„eiii,min cau be significantly lower than the ini- 
tial value of Ti in the outer belt, particularly after multiple 
scatterings backwards and forwards. Thus, this increases the 
range of potential orbits of scattered particles. 

This can be illustrated using an example system. Con- 
sider a particle scattered by the outer planet, with Ti = 2.99. 
The inner planet is placed arbitrarily at 02 = 0.7ai. The 
minimum pericentre for the particle after the particle is scat- 
tered by both planets, shown in Fig. |4l is qmin ~ 0.43ai. If 



5 MULTI-PLANET SYSTEMS 

All of the calculations discussed so far can be easily applied 
to planetary systems with many planets. The procedure dis- 
cussed in Sec. 14.11 can be repeated many times, to deter- 
mine the full range of orbital constraints and values for the 
Tisserand parameter after scattering by each planet. This 
analysis places useful constraints on the planets with which 
particles can interact, the planets that can eject particles 
and the minimum pericentre to which the whole system can 
scatter particles. 

All of the dynamics is determined by the initial value of 
the Tisserand parameter with respect to the outer planet, 
Ti, the outer planet's semi-major axis, ai and the ratio of the 
planets' semi-major axes to one another, Scaling the 

system, i.e. changing the semi-major axes, ai, whilst keeping 
their ratios, constant, will not affect the dynamics (val- 
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Figure 5. The variation in the minimum pericentre to which 
test particles can be scattered to by a system of five planets. 
The ratio of the planets' semi-major axes {a) is constant and is 
given as a ratio on the bottom axis and in terms of separation 
in Hill's radii, for five lOM^ planets, on the top axis. The initial 
value of the Tisserand parameter with respect to the outer planet 
is varied between 2.8 and 3.0. The shaded region illustrates the 
"unconstrained" regime for particles with T = 2.96, whilst the 
region to its left is the "non-interacting" and the region to its 
right is the "constrained" regime (see discussion in text). 

ues of Ti) and merely scales the minimum pericentre, qmin, 
with ai . In the next section, we discuss these constraints in 
terms of an example planetary system. 

5.1 A hypothetical 5-planet system with constant 
ratio of planets' semi-major axes 

We apply these calculations to a system of 5 planets, sep- 
arated by a constant ratio of adjacent planets' semi-major 
axis (-^^ = a). This corresponds to a constant number 
of Hill's radii for equal mass planets. Our results are inde- 
pendent of the planet masses. We fix the inner planet at 
as = aiN and calculate the semi-major axes of the other 
planets accordingly for a range of values for a. 

The minimum pericentre to which this system can scat- 
ter particles, shown in Fig. [S] as a function of a, is calcu- 
lated by repeatedly determining the minimum value of the 
Tisserand parameter for each planet. For the ith planet this 

occurs at g = gmin(Ti + l,min) (Eq. El, e = e;im,(ri+l,m,i„) 

(Eq. [U and cos / = 1. 

In this plot scattered particles exhibit three types of be- 
haviour. For simplicity we label the three types of behaviour 
as "non-interacting", "constrained" and "unconstrained". 
This refers to the constraints on the orbits of scattered par- 
ticles. In the "non-interacting" regime, the planets are so 
widely separated (small a) that particles cannot be scat- 
tered all the way along the chain of planets. The minimum 
pericentre to which one of the planets can scatter particles 
is outside of the next innermost planet's orbit. Hence the 
particles are restricted to the region surrounding the outer 
planet (s). 

In the "constrained" regime, the planets are so close to- 
gether (large a) that particles can be scattered between all 
planets in the system. If they are only scattered once along 



the chain of planets, the Tisserand parameter cannot vary 
significantly from its original value and there will be a non- 
zero minimum pericentre to which particles can be scattered. 
For such closely separated planets, it may no longer be valid 
to treat the scattering as a series of three body problems 
and the probability that particles are passed backwards and 
forwards between planets increases. This and the stability 
of planets so close together question whether particles scat- 
tered in any planetary system actually exhibit behaviour 
reminiscent of this "constrained" regime. 

As the separation of the planets is increased, the min- 
imum possible value of the Tisserand parameter for each 
planet decreases and hence the minimum pericentre for the 
whole system decreases. Eventually the separation is large 
enough that the Tisserand parameter falls below 2 and all 
constraints on the minimum pericentre are removed. This 
forms the third, "unconstrained" regime, where there are 
few constraints on the orbital parameters of scattered par- 
ticles. 

In Fig. [6] the constraints on the eccentricities and incli- 
nations of particles in the 3 regimes are shown. As particles 
are scattered by each planet, from the outermost (1) to the 
innermost (5), there will be a range of possible Tisserand pa- 
rameter values, between Ti,,„i»i (Eg. fT^ and Ti^rnax (Eg.fTl)) 
and hence a range of possible orbital parameters, given by 
Eqs. 1171 although of course the approximations used to cal- 
culate these may mean that they are not always strictly 
applicable, see discussion in ijTl It is the maximum incli- 
nation and eccentricity that are important on this figure, 
although of course the orbits of scattered particles will be 
distributed between the minimum and maximum values, in 
a manner not determined by this analysis. The plot shows 
that, for this example with Ti = 2.96, almost all planets can 
eject particles (e > 1) and that the scale height of the disc 
(inclinations of scattered particles) increases with decreas- 
ing distance to the star, as the constraints on the orbits of 
scattered particles decrease with each successive scattering 
event. It is clearly seen, as anticipated, that the constraints 
of orbits in the "constrained" regime are much tighter than 
those in the "unconstrained" regime. 

Although very few real planetary systems have planets 
separated by a constant ratio of their semi-major axes, it 
may be possible to similarly classify the behaviour of scat- 
tered particles within the three regimes and thus usefully 
better understand the future fate of scattered particles. 

5.2 Hypothetical multi-planet system separated 
by 107?// 

For real planetary systems the planets cannot be arbitrar- 
ily close togethe r as d ynamical instabilities are important. 
[chambers et al.l (| 19961 ) find that planets must be separated 
by at least 107?// to be stable. On Fig. (5] the separation 
of the planets is shown in terms of Hill's radii on the top 
axis, for a system of equal mass lOM^ planets. This shows 
that for the lOM^ planets considered, if they are separated 
by IQRh, then the behaviour of particles is unconstrained 
{qmin — >■ 0). Only very low mass (< lOM®) systems may be 
dynamically stable (separated by more than 107?//) and have 
limits on the scattering of particles, such that the particle's 
behaviour is in the "constrained" regime. 

Such low mass systems are, however, unlikely to only 
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Figure 6. Constraints on tiic eccentricities and inclinations 
(Eq. I17|l of particles scattered by a system of five planets with 
constant ratio of the planets' semi-major axes, a, and an initial 
value of the Tisserand parameter in the outer belt of Ti = 2.96. 
Particles are scattered from the belt, outside of planet 1, to the 
innermost planet, 5. Three planet separations are considered, 
corresponding to the three regimes (see discussion in the text); 
"non-interacting", a = 0.2, "unconstrained", a = 0.6 and "con- 
strained", a = 0.9. The dashed regions correspond to the param- 
eters of particles that can interact with the next interior planet 
(Eq. [9]l . The particles with high eccentricity were scattered out- 
wards and therefore are not on orbits that cross the inner planet's 
orbit. 



contain 5 planets. One possible outcome of planet forma- 
tion, is a chain of low mass planet embryos and an outer 
disc of planetesimals. Consider the example of such a disc in 
the position of the Solar System's Kuiper belt and a chain of 
interior, equal mass planets, between 1 and 30AU. If plan- 
ets Renerally form on orbi ts as tightly packed as possible 
l|Barnes fc Raym ond 2004; Raymond et al.ll2009t ). then their 
separation will be ~ 10i?H. We investigate the dynamics 
in such a system by varying the planet mass and thus the 
number of planets that fit between 1 and 30AU. This is 
equivalent to varying a. The results are shown in Fig. [T] 
The behaviour is identical to the five planet system in the 
"non-interacting" and "unconstrained" regimes, however the 
"constrained" regime no longer exists. This is because the 
interior planets further increase the parameter space avail- 
able to scattered particles. 



2.80 
2.88 
2.93 
2.96 
2.99 




0.001 



0.010 



0.100 1.000 10.000 100.000 1000.000 



Figure 7. The same as Fig. \5\ but for tightly packed planetary 
systems, with equal mass planets separated by lOii^f. The mass 
of the planets is shown on the bottom axis, whilst the top axis 
shows a. As many planets as fit between 1 and 30 AU are included, 
hence the minimum pericentrc is no longer finite for large a. 



6 APPLICATIONS TO REAL SYSTEMS 
6.1 Solar System 

This analysis can be applied to the planetary system that we 
understand best, our Solar System. There are three possi- 
ble sources of scattered bodies; the asteroid belt, the Kuiper 
belt and the Oort cloud. Most of the discussion so far has 
applied to the scattering of particles from a Kuiper-like belt, 
however, very similar processes occur in the asteroid belt. As 
discussed in iJJl we anticipate that the Tisserand parameter 
for objects scattered out of the Kuiper belt is close to 3. This 
should also apply to asteroids scattered from the main belt 
by Mars. The main difference between scattered asteroids 
and scattered Kuiper belt objects will be in the distribu- 
tion of the Tisserand parameter, TMars or Tjvep. This work 
does not determine these distributions, however, we specu- 
late that the distribution of TMars may be spread to lower 
values than that of Tncp- Jupiter is a strong perturber and 
may be able to alter the orbital parameters of an asteroid 
significantly in a single encounter. For Oort cloud comets 
scattered by planets, on the other hand, the range of val- 
ues of the Tisserand parameter is even larger, potentially 
including many comets with Tp < 2. This means that al- 
though this analysis is most usefully applied to scattered 
Kuiper-belt objects, it can also be applied to scattered as- 
teroids, but it cannot place many limitations on the orbits 
of scattered Oort cloud comets. 

Firstly considering objects scattered out of the Kuiper 
belt. Useful constraints can be made on their orbits in the 
outer planet region using this analysis; e.g. particle incli- 
nations are constrained to be below a maximum value, for 
example 80° , for scattered Kuiper belt objects with Tnbv > 
2.96, consistent with observations of Centaurs (|Gulbis et al.l 
It can also be inferred that the Solar System's outer 
planets are well placed for scattering particles between them. 
If Tjvep 2.982 then particles can be scattered, directly, 
all the way along the chain of planets to Jupiter and Ta- 
ble [2] shows that using Eq. [19] the planets are placed close 
to optimally for scattering particles as far inwards as pos- 
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Planet 


Semi-major 


axis (AU) 




Observed 


Optimum 


Neptune 


30.1 




Uranus 


19.2 


21.1 


Saturn 


9.58 


10.5 


Jupiter 


5.20 


3.3 


Mars 


1.52 




Earth 


1.00 


1.18 


Venus 


0.72 


0.74 


Mercury 


0.39 


0.33 



Table 2. The semi-major axes of solar system planets, compared 
the optimum semi-major axes in terms of scattering particles in- 
wards. These were calculated using fEg. I19| l for objects scattered 
from the Kuiper belt by Neptune with Tjvep = 2.98 and sepa- 
rately for objects scattered from the asteroid belt by Mars, with 

TMars = 2.98. 



sible. The three regimes presented in gSTT] can be applied 
to the Solar System to show that the majority of scat- 
tered Kuiper belt objects exhibit behaviour consistent with 
the "unconstrained" regime. These are shown in Fig. [S] If 
Tjvep < 2.962, there are no constraints on the orbits of 
scattered particles inside of Jupiter. On the other hand, if 
2.962 < Tjvep < 2.982 the dynamics of scattered bodies are 
"constrained" and there is a minimum value for the pericen- 
tre of objects scattered by Jupiter, whilst if TiVep > 2.982 
the particles are "non-interacting" and cannot be scattered 
into the inner planetary system. 

Similar constraints can be made for asteroids scattered 
into the inner planetary region. If Tuara < 2.97 then scat- 
tered objects behave as if they were in the "unconstrained" 
regime and if Tuara = 2.98, Eq.[T2]can be used in a similar 
manner to determine that the terrestrial planets are close to 
optimally separated for scattering particles between them, 
see Table[21 The important question, therefore, is what frac- 
tion of scattered of asteroids have Tuars ~ 2.98. If TMars 
is lower then issues arise with some of the approximations 
used in this work. For example, particles may not be scat- 
tered directly along the chain of planets. For Tuars < 2.96, 
the minimum pericentre to which asteroids may be scattered 
by Mars is already inside of Venus's orbit, qmin < avenua 
(Eq. [31). Particles may interact directly with Venus before 
being scattered by Earth, or be scattered multiple times 
by both Venus and Earth. As described in H4.4I this would 
greatly increase the range of potential orbits for scattered 
particles. Also mean motion resonances and secular effects 
will play a greater role in altering the dynamics of scattered 
asteroids than of scattered Kuiper belt objects, raising fur- 
ther doubt as to the validity of this approach. It can nonethe- 
less be used cautiously to further investigate the dynamics 
of asteroids scattered into this region. 



6.2 Warm dust discs 

As discussed in the Introduction, there are many obser- 
vations of stars with warm dust belts, e.g. (|Gaidosl 1 19991 : 
iBeichman et all l2005l : ISong et al.l 120051 1. Many of the sys- 
tems with warm dust also have cold dust belts, amongst 
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Figure 8. The minimum pericentre (qmin) to which Kuiper 
belt objects can be scattered to by Neptune and the outer Solar 
System planets (solid line) and similarly for asteroids scattered 
by the terrestrial planets (dashed line), as a function of the ini- 
tial value of the Tisserand parameter with respect to Neptune or 
Mars, respectively. The shaded area illustrates the "constrained" 
region, for scattered Kuiper belt objects with diagonal shading 
in green and for scattered asteroids with vertical shading in red. 
The "unconstrained" region lies to the left of the plot and the 
"non-interacting" regime to the right. 



others, r? Corvi ([Smith et al.ll2009l:IWvatt et al.ll2005l). /3Leo 
jChurcher et allboill ) and e Eri (|Backman et al.ll2009l ). The 

analysis presented here can be used to consider the scatter- 
ing of particles from an outer belt inwards, as a potential ex- 
planation for the observed warm belts. Our main conclusion 
is that the architecture of a planetary system determines 
whether or not material can be scattered to the position of 
the observed belt. We, therefore, speculate that the diversity 
of planetary system architectures could result in the diver- 
sity of observed systems, both in terms of disc radii and the 
ratio of the flux from the outer to the inner belt. Although 
this analysis does not determine what fraction of the scat- 
tered particles end up in the position of the observed disc, it 
does show that some planetary systems cannot scatter par- 
ticles onto the required orbits and illustrates that when the 
distribution of scattered particles is determined, tight con- 
straints will be placed on the architecture of the planetary 
system required. 

Consider the example system of 77 Corvi, with cold and 
warm du st. The inner belt i s resolved and lies between 0.16- 
2.98 AU (|Smith et allbOOgl) whilst the outer dust is at 150 
± 20 AU (|Wvatt et all 120051 ). Ahhough there are no plan- 
ets detected in this system, it seems probable that there is 
a planet clos e to the inner edge of the cold, outer belt, that 
truncates it (|Wvatt et al.ll2005l ). We, therefore, consider a 
planet at lOOAU. If the Tisserand parameter with respect 
to this planet is Ti = 2.96, then this planet alone could po- 
tentially scatter particles in as far as 47AU (Eq. [3]). In order 
for particles to be scattered inwards to the location of the 
warm belt, Qmin < 3AU, at least three planets are required. 
The optimum position for these planets is 58 and 23 AU, 
with the outer planet at 100 AU (Eq. I19[). The orbits cannot 
vary significantly from these values if the minimum pericen- 
tre is to remain less than 3 AU. For example, if the planets 
were positioned at 100, 80 and 60 AU, particles could only 
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be scattered in as far as 6AU and thus the warm dust beh, 
if it formed, would be at larger radii. Alternatively, there 
could be more than 3 planets, the initial value of the Tis- 
serand parameter could be less than 2.96 or particles could 
be scattered multiple times backwards and forwards between 
the planets, as discussed in ^4.41 

6.3 Metal polluted white dwarfs and white dwarfs 
with close-in circumstellar discs 

Evidence of evolved planetary systems and scattering of 
planetary material is found in the observations of metal pol- 
luted white dwarfs (jZuckerman et aP l2003l : iKoester et al.l 
I2OO5I) and white dwarfs with close-in circumstellar discs 
( Farihi et al]|2009l ). In order to explain these observations 
with planetary material, comets or asteroids must be scat- 
tered onto star-grazing orbits and tidally disrupted. The 
analysis presented in this work can be used to determine 
the feasibility of this explanation. 

Planets are required to scatter comets or asteroids close 
enough to the star. There are three potential reservoirs in an 
evolved planetary system, a Kuiper belt analogue, an Oort 
cloud analogue and if it survives an asteroid belt analogue. 
This analysis shows that it is possible for particles from all 
three reservoirs to be scattered onto star-grazing orbits, but 
that this ability depends strongly on the planets' orbits and 
the initial value of the Tisserand parameter. The lower the 
initial value of the Tisserand parameter, the more likely that 
particles can be scattered sufficiently close to the star (the 
lower Qmin Fig. [2]). Hence, the majority of particles from an 
Oort cloud analogue can be scattered onto star-grazing or- 
bits, whilst for a Kuiper or asteroid belt analogue this ability 
is strongly dependent on the initial value of the Tisserand 
parameter and the planets' orbits. Asteroid belt analogues 
have the advantage of lower initial values for the Tisserand 
parameter, but the disadvantage that there may be fewer 
surviving interior planets and the asteroid belt itself may 
not survive until the white dwarf phase. 

There are a large number of observations of Kuiper belt 
analogues around main sequence stars (jWvatt 200i ) and 
models find that such systems survive the star's evolution 
(IBonsor fc Wvatt|[2O10l ). Such belts have been su ggested as 
the source of the metal pollution e.g. I Jural (|2008l ). although 
there is little evidence that they are capable of scattering 
particles sufficiently close to the star. Here, we show that it 
is possible for some planetary systems to scatter particles 
from an outer belt onto star-grazing orbits, but that there 
are tight constraints on the planets' orbits and the initial 
value of the Tisserand parameter in the outer belt. 

One potential hindrance in the ability of an evolved 
planetary system to scatter particles onto star-grazing or- 
bits is the absence of i nner planets due to the star's evolu- 
tion. IViUaver fc Lhdol (|2007l i find that white dwarfs should 
not possess planets within 15AU due to a combination of 
the increased stellar radius, tidal forces and stellar mass 
loss. In order for a planet at at = 15AU to scatter parti- 
cles onto star-grazing orbits {qmin < Rq), particles must 
have values of the Tisserand parameter less than 2.05 when 
they interact with the planet (Eq. [3|. Only particles from 
an evolved Oort cloud might have sufficiently low values 
of Tisserand parameter without interacting with further 
planets. Therefore, using repetition of the technique de- 



scribed in i]4.1l if particles originate in an outer belt with 
Ti > 2.97, then at least 4 planets are required to scatter par- 
ticles onto star-grazing orbits, whilst for 2.89 < Ti < 2.97 
only 3 are required. Another potential hindrance is the in- 
stability of many planet ary systems after stellar m ass loss 
on the giant branches (|Debes fc SigurdssonI |2002| ) , if, for 
example, planets are ejected. Examples of real planetary 
systems that could scatter particles onto star-grazing or- 
bits from a Kuiper-like belt include our Solar System (if 
TNep < 2.96) and HR 8799 with planets at 14.5, 24, 38 and 
68 AU jMarois et al.ll2008l . lioiol '). if Ti < 2.95 in the outer 
belt. 

This analysis crucially shows that it is possible to scat- 
ter comets or asteroids onto star-grazing orbits and places 
limits on the architecture of a planetary system that can do 
this, but it does not inform us about the probability of a 
given planetary system to scatter planetesimals onto star- 
grazing orbits. Oort cloud analogues only require a single 
planet to scatter material onto star-grazing orbits, whilst 
constraints are placed on the orbits of planets and the initial 
value of the Tisserand parameter required to scatter mate- 
rial inwards from a Kuiper or asteroid belt analogues. Thus, 
this analysis shows that material from an evolved Kuiper 
belt is a potential origin of the metal pollution in white 
dwarfs, although fewer constraints exist on the ability of 
an evolved Oort cloud to scatter comets onto star-grazing 
orbits. This provides important evidence in support of the 
planetary origin for the white dwarf observations. 



7 DISCUSSION OF LIMITATIONS 

The purpose of this work is to present a simple analytical 
tool that can be applied to many planetary systems. It de- 
termines the potential orbital parameters of scattered par- 
ticles, based on the initial value of the Tisserand parameter 
and the planets' orbits. It does not claim to determine the 
probability for any particle to be scattered onto a given or- 
bit, nor the expected distribution of scattered particles. In 
order to retain this simplicity it was necessary to make sev- 
eral assumptions that strictly limit the applicability of this 
analysis. These are discussed below. 

We anticipate that the behaviour described in this work 
will be useful in interpreting N-body simulations of small 
body scattering in planetary systems. In such simulations 
the initial value of the Tisserand parameter for individual 
particles will be well constrained, and this means that it 
should be possible to describe their subsequent evolution us- 
ing the analysis presented here. However, this presupposes 
that our simplifying assumptions are valid, and it will be 
important to use N-body simulations to test this. The main 
assumptions regarding the dynamics that we see are sum- 
marised here. 

One of the biggest limitations in the analysis presented 
here, is the dependence on the initial value of the Tisserand 
parameter in the outer belt. This is in general an unknown 
quantity, although good approximations can be made to its 
value, as discussed in ij2] and its value is determined in N- 
body simulations. 

Here, the scattering of particles by a chain of planets is 
considered as a series of three-body problems. This should 
be broadly true, although some particles may be affected by 
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Pericentre 

Figure 9. The extension of the orbital parameter space available 
to scattered particles with Tp = 2.98. This is analogous to the 
plots in the top row of Fig. [T] except that the particles are allowed 
to interact with the planet within a zone of size A ~ Rh^ for 
Jupiter mass (green) and Neptune mass (red) planet. For the 
Neptune mass planet the orbital parameter space available to 
scattered particles is not altered significantly from that shown in 
Fig. [T] whereas for the Jupiter mass planet the strict limit on the 
minimum pericentre is removed. 



secular or resonant perturbations, or interact with a planet 
other than the one dominating their dynamics during that 
period. This could alter the value of the Tisserand param- 
eter. Particles may also be passed backwards and forwards 
along the chain of planets, as discussed in ii4.4l 

Another very important limiting assumption is that 
particles only interact with the planet when their orbits ex- 
actly cross the planet's orbit, whereas in reality there will be 
a small zone around the planet within which a particle may 
be scattered by the planet. The size of such a zone depends 
on the planet's mass and therefore its inclusion into the cal- 
culations would introduce a mass dependence to the analysis 
presented here. For the low mass planets all of the analysis 
presented here should be valid, whilst for higher mass plan- 
ets care must be taken when rigorously applying some of the 
derived limits. This is illustrated in Fig. |5] which shows the 
extension to the orbital parameter space for particles that 
can interact with the planet in a zone of size — . This ex- 

^ ap 

tends the range of potential orbital parameters of scattered 
particles, as orbits with Q > ap(l — A) and q < ap{l + A) 
can interact with the planet, rather than just Q > ap and 
q < Up. This orbital parameter space is shown in Fig. [9] for 
Jupiter and Neptune, with A = Rh and Tp — 2.98. Jupiter 
is massive enough that there is no longer a limit on the min- 
imum pericentre to which particles can be scattered and the 
available parameter space is increased by a small, but sig- 
nificant, amount. Neptune, on the other hand, does not sig- 
nificantly extend the parameter space available to scattered 
particles and therefore all of the analysis presented in this 
paper will apply. For Jupiter mass planets the broad results 
presented here can still be applied, however, care should be 
taken when rigorously applying the derived limits. 

Strictly the conservation of the Tisserand parameter, 
and therefore this analysis, should only be applied to sys- 
tems with co-planar planets on circular orbits, i.e. within 



the context of the circular restricted three body problem. 
It is, however, found that even when these assumptions are 
rel axed, the analysis st i ll app lies approximately, for exam- 
ple lMurrav fc DermottI (|l999i ) found only a small change in 
the Tisserand parameter when they consider Jupiter's eccen- 
tricity. Caution should, however, be exerted when applying 
this analysis to some of the detected exoplanets with large 
eccentricities (and relative inclinations). 



8 CONCLUSIONS 

The purpose of this work was to describes simply and an- 
alytically the scattering of particles in any planetary sys- 
tem. This analysis constrains the outcomes of scattering 
events, based on the conservation of the Tisserand param- 
eter (Eq. in a manner that is very useful for analysing 
the structure of many planetary systems where the scatter- 
ing of small bodies by planets is important. This analysis 
can be used to better understand behaviour seen in N-body 
simulations. 

We consider here the application to planetary systems 
where small bodies are scattered from an outer belt by inte- 
rior planets. The analysis could, however, easily be reformu- 
lated to consider scattering by planets exterior to the belt. 
We assume that the scattering process can be approximated 
by a series of three-body problems, during each of which the 
Tisserand parameter with respect to the relevant planet is 
conserved. This constrains the possible orbits of scattered 
particles, based solely on the initial value of the Tisserand 
parameter and the orbits of the planets, with the assumption 
that particles are passed directly along the chain of planets 
and that particles only interact with the planet when their 
orbits directly cross the planet's orbit. In this case there is 
no dependence on the planet's mass and it is only the ratio of 
the planets' semi-major axes that are important. A depen- 
dence on planet mass would be introduced if the interactions 
of particles with the planet within a zone around the planet 
were included. This analysis places an important limit on 
how far in particles can be scattered {qmin from Eq. [3]) and 
determines which planets the particles can interact with, 
which can eject them and the potential height of the disc, 
based on the maximum particle inclinations (Eq. I17p . We 
consider the full range of possible orbits of scattered parti- 
cles, rather than their distribution. 

In this work we consider the application of this analy- 
sis to our Solar System, main sequence stars with both cold 
and warm dust belts and metal polluted white dwarfs. In 
the Solar System, this analysis describes the scattering of 
Kuiper belt objects by Neptune to become Centaurs and 
Jupiter Family comets, as well as asteroids by Mars and the 
terrestial planets. We show that the Solar System planets 
are close to optimally separated for scattering particles be- 
tween them. One explanation for main sequence stars with 
warm dust belts that cannot have survived for the age of 
the system in their current positions is that material was 
scattered inwards from an outer belt. If this is the case, this 
analysis shows that certain architectures for the planetary 
system could not produce the observations. Given the strong 
dependence of the scattering process on planetary system ar- 
chitecture, we speculate that the diversity of such systems is 
a reflection of the variety of planetary system architectures. 
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Observations of metal polluted white dwarfs and white 
dwarfs with circumstellar discs have been associated with 
material scattered inwards from an outer evolved planetary 
system. Such material can be scattered sufficiently close to 
the star from all Oort cloud analogues and, given certain 
constraints on the system architecture, from some Kuiper 
belt analogues. This strengthens the case for a planetary 
origin to these observations, although this analysis does not 
comment on the probability for particles to be scattered onto 
such orbits. 

In summary, the analytical tool presented here can aid 
our understanding and place useful constraints on the scat- 
tering of small bodies in a wide range of planetary systems. 
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